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Abstract 

Exploiting spectral properties of symmetric banded Toeplitz matrices, we de- 
scribe simple sufficient conditions for positivity of a trigonometric polynomial for- 
mulated as linear matrix inequalities (LMI) in the coefficients. As an application 
of these results, we derive a hierarchy of convex LMI inner approximations (affine 
sections of the cone of positive definite matrices of size m) of the nonconvex set of 
Schur stable polynomials of given degree n < m. It is shown that when m tends 
to infinity the hierarchy converges to a lifted LMI approximation (projection of an 
LMI set defined in a lifted space of dimension quadratic in n) already studied in the 
technical literature. 
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1 Introduction 

Linear system stability can be formulated algebraically in the space of coefficients of the 
characteristic polynomial. The region of stability is generally nonconvex in this space, and 
this is a major obstacle when solving fixed-order or robust controller design problems. 
In the case of discrete-time linear systems, the region of stability is a bounded open 
set whose boundary consists of (flat) hyperplanes and nonconvex (negatively curved) 
algebraic varieties. Recent results on real algebraic geometry and generalized problems 
of moments can be used to build up a hierarchy of convex linear matrix inequality (LMI) 
outer approximations of the region of stability, with asymptotic convergence to its convex 
hull, see e.g. |7j for a software implementation and examples. It is generally more difficult 
to construct convex LMI inner approximations, see [5] for a survey. Strict positive realness 
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of rational transfer functions and its connection with polynomial positivity conditions are 
used in [BJ to generate inner approximations which are lifted LMI sets. For polynomials 
of degree n, they are projections onto coefficient space M. n of an LMI set living in a lifted 

space R 2 . The LMI set is built around a particular point, the central polynomial, 
whose relevance in robust control design is explained in [BJ. These lifted LMI regions 
are also used in signal processing, see e.g. [31 Section 7.3]. They can be derived in a 
state-space setting with the Kalman-Yakubovich-Popov lemma @]. 

Whereas lifted LMIs are a powerful modeling paradigm (it is currently conjectured that 
every convex semialgebraic set is a lifted LMI set), the introduction of a large number 
of lifting variables can be seen as a drawback. It is therefore relevant to build convex 
LMI inner approximations of the nonconvex stability region without liftings, namely as 
affine sections of the cone of positive semidefinite matrices. This is the objective of 
this paper. We use results of functional analysis on sequences of eigenvalues of Toeplitz 
matrices to derive sufficient LMI conditions for positivity of trigonometric polynomials, 
and we apply these results to construct a hierarchy of m-by-m LMI inner approximations 
of the nonconvex stability domain. Moreover we prove that when m tends to infinity, the 
hierarchy converges asymptotically to the lifted LMI approximation of [BJ. 



2 Trigonometric polynomials and Toeplitz matrices 



Let pk,k = 0,l,2,...,n denote real numbers, and define the trigonometric polynomial 

z = e ie h> p(9) =p + Pl (z + z~ l ) + p 2 (z 2 + z~ 2 ) H h p n (z n + z~ n ) 

= Po + 2pi cos 6 + 2p 2 cos 26 + ■ ■ ■ + 2p n cos n6 

of degree n mapping the unit disk of the complex plane onto the real axis. 

For a given integer m > n, define the column vector v m (z) = [1 z z 2 ■ ■ ■ z m ~ l Y 
represent polynomial p as a quadratic form 



and 



where 



p(9) = -vl(e- ie )P m v m (e ie ) 



Po — x Vi 

m m 



PO 



is an m-by-m symmetric banded Toeplitz matrix. 
Define 

Po Pi P2 
Pi Po Pi 

R m = P2 Pi PO 



Po 



(1) 



Po 



(2) 



2 



as the m-by-m moment matrix of p, so named for 

»2tt 



Pk 



2tt 



p(8)e- lk "d8 



o 



is the k-th moment, or Fourier coefficient, of polynomial p. Note that R m has the same 
banded symmetric Toeplitz structure as P m . 

Connections between the spectrum of matrix R m and the values taken by polynomial p 
on the unit circle have been studied extensively. In the sequel, A m i n denotes the minimum 
eigenvalue of a symmetric matrix. 



Theorem 2.1 



lim A min (P m ) = mmp(6). 



Proof: It is a corollary of Gabor Szego's fundamental eigenvalue distribution theorem, 
see e.g. [3 Corollary 4.2]. □ 

In this section we aim at establishing a similar spectral property linking matrix P m and 
polynomial p. First let us state a few instrumental results. 



Lemma 2.1 For all 9 it holds A m i n (P m ) < p(9) and as a consequence 

lim sup A min (P m ) < min p(6). 



(3) 



Proof: From relation ([T]) and the identity v^ie )v m (e ) = m, it follows that 



m 



i£(e-")P m i; m ( C ") 



p{9) 



and hence A m i n (P m ) < p{6). When ra-^oowe obtain the desired result. 
Lemma 2.2 



0{m 



Proof: Consider 



Rr, 



1 



m-lP 1 



-Pi 







1 



m— 1 



rn—1 

—^P2 —^—rPi 



Pi 



and hence for the Froebenius norm 



IP - R I 



E 



m — k 



J (m — k) 



- 1 

iVk = ^ k 

k=l 



Pk- 



We are now ready to state our main result. 



(4) 
□ 



□ 
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Theorem 2.2 



lim A min (P m ) = minp(0). 

m— >+oo 8 



Proof: Let v be an eigenvector of P m such that v T v = 1 and P m v = A m i n (P m )f . From 
the equality 

v T P m v = v T (P m - Rm)v + v T R m v, 
we obtain with the help of Lemma 12.21 the following inequality 

Amin {P m ) > 0(m 2 ) + A min (P m ). 

Taking the limit, we obtain 

liminf A min (P m ) > lim A min (P m ). 

m— >+oo m— >+oo 

Using Lemma 12.11 and Theorem I2.1[ we can see that 



and hence 



liminf A min (P m ) > lim \ mm (R m ) = mm p(6) 

m— >+oo m— >+oo 8 



liminf A m i n (P m ) > mmp(9) > lim sup A m i n (P„ 

m-5-+oo 8 m->+oo 



If x m is a real sequence then it is well-known that if lim inf m _s. +00 x m 
then the sequence x m converges to \im m ^ +00 x m = lim inf m ^ +OC) x m 
and this completes the proof. 



lim sup m _ > _|_ 00 x m , 



m— >+oo ^mi 
□ 



lim sup 



Corollary 2.1 Assume that polynomial p is positive. Then, there exists a sufficiently 
large integer m such that for m > m , the Toeplitz matrix P m is positive definite. 



Proof: Use Theorem 12.21 



□ 



Remark 2.1 Note that when p is positive, matrices P m are not necessarily positive def- 
inite if m is not large enough. As a simple example consider the positive polynomial 
p{9) 



2 + 2cos#+ |cos26». We have 



9 3 12 

Z 2 5 

3 o 3 

2 Z 2 

12 3 9 

5 2 Z 



which is not positive definite, since A m i n (P3) = — |. Also, the next Toeplitz matrix 



P, 



2 ! ! o 

4 9 4 8 

3 Z 3 5 
8 4 9 4 

5 3 Z 3 

o ! ! 2 



is not positive definite either, since A m i n (P4) 



_8_ 

13 



2V509 
15 



—0.3415. However, one can 



check that when m > m = 30, matrices P m are indeed positive definite. 



4 



3 LMI inner approximations of stability domain 



Consider a monic polynomial 

d(z) = do + diZ-l h d n „iz n ~ l + z n 

of degree n, with coefficient vector d G lR n and let us define the set 

S = {deR n : d(z) stable} 

where stability is meant in the discrete-time, or Schur sense, i.e. all the roots of d(z) belong 
to the open unit disk. Many control problems (e.g. fixed-order or robust controller design) 
can be formulated as linear programming problems in S. Unfortunately S is nonconvex 
when n > 2, which renders controller design difficult in general. It can therefore be 
relevant to describe convex inner approximations of S, in particular by exploiting the 
modeling flexibility of linear matrix inequalities (LMIs), see [5] and references therein. 

An approach consists in choosing a monic polynomial 

c(z) = C + CiZ H h Cn^iZ^ 1 + Z n 

which is stable. Once c is given, we define the trigonometric polynomial 

z = e ie ^ p c ' d (9) = c(z- l )d(z) + c(z)d(z- 1 ) 



2 Cjdf. cos 19 



1=0 j,k=0 

\j-k\=i 



and the set 



V c = {d G R n : p c ' d {6) > V 9 G R}. 
Lemma 3.1 Let c(z) be a given stable polynomial. Then V c C S. 

Proof: A geometric proof is as follows. Since polynomial c(z) is Schur stable, when z = e ld 
varies along the unit circle, complex number c(e %e ) has a net increase of argument of 2mr, 
or equivalently the plot of c(e ld ) encircles the origin n times, see e.g. [21 Section 1.3.3] or 
use Cauchy's argument principle. Notice that the real number p c,d (9) = c(e~ ld )d(e ld ) + 
c(e z6 )d(e~ td ) is equal to 2\c(e zd )d(e td )\ cos(c(e* e> ), d(e z6 )) where the last term is the cosine 
of the oriented angle between vectors c{e td ) and d(e ) in the complex plane. Therefore 
p c,d (9) positive implies that the cosine is positive and hence that the angle between c(e t6 ) 
and d(e l9 ) is less than | in absolute value for any given value of 9. This means that 
complex number d(e %e ) also encircles the origin n times when 9 range from to 2ir, and 
hence that polynomial d(z) is Schur stable. □ 

Let P£f be the symmetric banded Toeplitz matrix corresponding to polynomial p c,d , built 
as in (j2J), and define the set 

V c m = {d G R n : P c m d y 0} 

where >~ means positive definite. Note that symmetric matrix P^ d depends affinely on 
d, so that T"^ is a convex LMI set. 
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Theorem 3.1 Let c(z) be a given stable polynomial of degree n, and let m > n. Then 
V°cS. 



m 



Proof: Since mp c ' d {6) = i>^(e %e )P^ b d v m {e ld ) 1 positive definiteness of matrix P^f implies 
positivity of polynomial p c,d (8). Then use Lemma [3. II □ 

Set is therefore a valid convex inner approximation of the nonconvex stability region 
S. Its geometry depends only on the choice of a stable polynomial c(z). 

Theorem 3.2 Let c(z) be a given stable polynomial. Then V c = \im. m ^ +ao l :> ^ n . 

Proof: Use Theorem O □ 

Finally we make the connection with the results in [6\. Recall that a discrete-time rational 
function is strictly positive real (SPR) whenever its real part is strictly positive when 
evaluated along the unit circle. 

Theorem 3.3 

V c = {d G W n : 44 SPR}. 

c{z) 

Proof: Since c(z) is stable, the SPR inequality 

d{e ie ) _ 1 ( d(e*) d{e- i0 ) \ _ c{e- ie )d{e ld ) + c{e ie )d{e~ ld ) 

is equivalent to positivity of trigonometric polynomial p c,d (8). □ 

Polynomial c(z) is referred to as a central polynomial in [6] since set V c is built around 
c(z) in the coefficient space. Note however that there is no guarantee that c(z) belongs 
to if m is not large enough, see Remark 12.11 



4 Example 



4.1 Second-order polynomials 



We consider second-order polynomials for which the exact stability region is a triangle 
with vertices (z + l) 2 , (z — l)(z — 1) and (z — l) 2 [H Example 11.13]. 

Choosing c(z) = z 2 , we have p c,d {6) = 2 + 2d\cos6 + 2docos29. The first LMI inner 
approximation is 



V^ d = {(do,d 1 ) : P, 



c,d 



2 |dj 3d, 



Mi 



2 5 

Z 2 







3d |d a 2 



^0} 



and it is represented on the left of Figure [T] within the stability triangle, as claimed by 
Theorem 13.11 
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Figure 1: 3-by-3 LMI set (shaded gray, left) and 4-by-4 LMI set (shaded gray, right) 
within second-order discrete-time stability region (triangle). 



The second LMI inner approximation is 



c,d 



{(do,di 



P 



c,d 



2 

2d 




\d\ 2d[ 



o 



-A 

2d 





^d\ 2do 

4 
3 



2 U x 



^0}, 



see the right of Figure [TJ 

On the left of Figure [2] we represent the LMI set 



c,d 



2 

5^0 








\d x 


7 f 











2 


|di 
2 


I do 
2 








Mi 







l d o 


U x 


\d 





Uo 


Idi 


2 










Ido 




2 













¥ 








Ido 

Ui 
2 



yo}. 



The boundary of this set is piecewise polynomial, defined by two algebraic plane curves 
whose irreducible defining polynomials -7200 + 50A0d + 3528dl + 4900^ - hlAhd d\ (a 
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Figure 2: 7-by-7 LMI set (shaded gray, left) and 50-by-50 LMI set (shaded gray, right) 
within second-order discrete-time stability region (triangle). 



cubic) and 6480000 + 4536000d - 9525600^ - 8820000d 2 



4445280dg 



7717500d d? + 



3111696^ 

ic.d 



4321800dgdf + 1500625df (a quartic) factor the determinant of the 7-by-7 
pencil Pj' a . See Figure |3] for a representation of this set and the algebraic components of 
its boundary. 



On the right of Figure |2] we represent the LMI set V§ which, according to Theorem 13.2 
is almost equal to the lifted LMI set 



V c ' d = {(cMO : 3(g ,gi,<72) : 



go gi d 
gi Q2 - go di- gi 
d d x - gi 2 - g 2 



^0}, 



the projection onto R 2 of an LMI living in M 5 , and which is the union of an ellipse and a 
triangle, as studied in [6]. 



4.2 Third order 

We consider third-order polynomials for which the exact stability region is delimited 
by a nonconvex hyperbolic parabolic embedded in a tetrahedron with vertices (z + l) 3 , 
(z + l) 2 (z - 1), (z + l){z - l) 2 and (z - l) 3 , see [H Example 11.14]. 
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Figure 3: 7-by-7 LMI set (shaded gray) and the algebraic components of its boundary 
(thick black lines). 



Choosing c(z) 



we have p c ' d (9) = 2 + 2d 2 cos 9 + 2d 1 cos 26 + 2d cos 39. The first LMI 



inner approximation is 



V c / = {(d ,d 1 ,d 2 ) : P L 



cA 



3^ 2 

2d 1 
Ad 



\d 2 2d 1 
2 Ui 



4 rJ 

3 2 

2d x 



\d 2 



4d 
2di 
\d 2 



^0} 



and it is represented on the left of Figure H] within the nonconvex stability region, as 
claimed by Theorem 13. 1[ On the right of Figure H] is represented the LMI set V^q which, 
according to Theorem 13.21 is almost equal to the lifted LMI set V c,d . 



5 Concluding Remarks 

We have used results on spectra of Toeplitz matrices to construct a hierarchy of convex 
inner approximations of the nonconvex set of stable polynomials, with potential applica- 
tions in fixed-order robust controller design. The main difference with respect to previous 
results is that the inner sets are defined by LMIs (affine sections of the cone of positive 
definite matrices) without the need to resort to projections and lifting variables. More- 
over, our LMI sets belong to a hierarchy converging asymptotically to a lifted LMI inner 
approximation described previously in [6]. 
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Figure 4: 4-by-4 LMI set (shaded gray, left) and 50-by-50 LMI set (shaded gray 
right) within third-order discrete-time stability region (delimited by a meshed hyperbolic 
parabolic embedded in a tetrahedron). 



It is likely that our results can be extended to deal with positive trigonometric polynomial 
matrices and block Toeplitz matrices, with potential applications in multi-input multi- 
output control systems. 

Sufficient conditions ensuring that a real polynomial is a sum-of-squares (and hence that 
it is positive) have been proposed in [8J, so it could be insightful to transpose these 
conditions to trigonometric polynomials and compare with our approach. Results in [8] 
are also valid for multivariate polynomials, and this may have applications in fixed-order 
or robust controller design for multi-dimensional systems. 
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